Optimal quantum tomography of permutationally invariant qubits 
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Mutually unbiased bases determine an optimal set of measurements to extract complete information about 
the quantum state of a system. However, quite often a priori information about the state exist, making some 
of the measurement bases superfluous. This is, for example, the case when multiqubit states belong to the 
permutationally invariant subspace. In this paper we derive the minimal sets of mutually unbiased bases needed 
to tomographically reconstruct such states. 
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I. INTRODUCTION 

The quantum state is a mathematical entity that encodes 
complete information about a system: once it is known, the 
probabilities of the outcomes of any measurement can be pre- 
dicted. It seems thus indisputable that ascertaining an un- 
known state accurately turns out to be of uttermost importance 
for modern quantum technologies. Broadly speaking, this is 
the scope of quantum tomography |! l] which, over the past 
years, has evolved from the initial theoretical |2| and experi- 
mental concepts y[ to a widely acknowledged and fairly stan- 
dard method used extensively for both discrete |4] and contin- 
uous JH variables. 

However, the tomographic task becomes harder as we ex- 
plore more intricate systems. For example, for the simple case 
of n qubits, 2 2 ' 1 — 1 real numbers are required for its complete 
characterization, while any von Neumann measurement gives 
only 2" — 1 independent data. Consequently, one will have to 
make at least 2" + 1 different such measurements before one 
can claim to know everything about an a priori unknown sys- 
tem. With such a scaling, it is clear that the methods rapidly 
become intractable for present state-of-the-art experiments. 

As a result, more sophisticated tomographical techniques 
are called for. New protocols try to exploit the idea that the 
scheme is explicitly optimized only for a particular kind of 
states [6]. In that perspective, we look here at the specific but 
not unimportant example of n qubits prepared in an arbitrary 
state that is, however, known to be invariant with respect to 
any qubit permutation. This may be due to, e.g., a permuta- 
tionally invariant preparation Hamiltonian. In this instance, 
the associated Hilbert space has dimension « + 1, and there- 
fore it should be possible to reconstruct such a state with only 
n + 2 von Neumann measurements. 

Permutationally invariant qubit states are employed in di- 
verse quantum information strategies 0. They are also opti- 
mal for quantum metrology 0] and play an important role in 
the characterization of locally non-interconvertible entangle- 
ment classes Through all this paper we take permutational 
invariance for granted; theoretical tests of this property (others 
than full tomography) have been put forward in [10s]. 

Recently, a number of suggestions have appeared for an ef- 
ficient generation of different entangled permutationally sym- 



metric qubit states ITTll . The tomography of such states has al- 
ready been discussed in Ref. lfl2tl . and a four-qubit experiment 
has been performed lfl3ll . However, in these proposals the 
measurements have been chosen as a set of informationally 
complete projectors. This may provide a certain experimental 
simplicity (e.g., for spin states it may be possible to simply 
use the orientation of a Stern-Gerlach apparatus to chose the 
projector), but is by no means an optimal strategy. 



The number of separate von Neumann measurements 
needed for a complete state determination is optimal when the 
bases in which those measurements are performed are mu- 
tually unbiased ifLHl (in the standard n qubit Pauli tomogra- 
phy, 3" different settings are needed, while in this optimal 
approach, 2" + 1 are enough). In fact, the notion of mutu- 
ally unbiased bases (MUBs) emerged in the seminal work of 
Schwinger |fl5}l and has turned into a cornerstone of quantum 
information, mainly due to the spotlight placed on them by the 
elegant work of Wootters and coworkers |fl6l| . MUBs are en- 
dowed with the property of being maximally incompatible, in 
the sense that a state giving precise results in one set (i.e., one 
of the basis states) produces maximally random results when 
measured in another basis in the MUB set. 



Another remarkable advantage of the MUB-based tomog- 
raphy is that each measured probability determines a single 
element of the density matrix so, in principle, there should be 
no need for a numerical data inversion to reconstruct the corre- 
sponding state. In practice, however, experimental noise and 
measurement imperfections may yield an unphysical density 
matrix, so a fitting procedure might still be needed. 



For all these compelling reasons, we think it is worthwhile 
to first prove that minimal complete sets of MUBs exist for the 
tomography of a permutationally invariant n-qubit state, and 
subsequently show how to construct them. This is precisely 
the goal of this paper. Of course, these advantages come with 
a price: those minimal MUB sets are, in general, comprised of 
entangled projectors, which renders their experimental imple- 
mentation more challenging than the product Pauli projectors. 
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II. MUTUALLY UNBIASED MEASUREMENTS FOR 
QUBITS 

A compact way of labeling «-qubit states consists in using 
the finite field F2" (the reader interested in mathematical de- 
tails is referred, e.g., to the excellent monograph by Lidl and 
Niederreiter lfl7ll ). For our purposes, this can be considered 
as a linear space spanned by an abstract basis {6\ , . . . , 0„}, so 
that given a field element v (henceforth, they will be denoted 
by Greek letters) the expansion v = YH=\ n i ®i (with «; G Z2) 
allows us the identification v 44> (ni,«2, • • • ,n„). 

Moreover, the basis can be chosen to be orthonormal with 
respect to the trace operation (the self-dual basis); that is, 
Iv(djdj) = 8ij, where tr(v) = v + v 2 + . . . + v 2 " ' and maps 
¥2" onto the base field Z?. In this way, to each qubit we as- 
sociate a particular element of the self-dual basis: /th qubit 
0, 

Let {|v)} be an orthonormal basis in the Hilbert space of 
the system, which is isomorphic to C 2 ". Operationally, the 
elements of the basis are labeled by powers of a primitive el- 
ement. These vectors are eigenvectors of the operators Z a 
belonging to the generalized Pauli group llHHl . whose basic 
generators are 

Z a = £(-l)«'<™>|v)<v|, Xp=£|v + /3)(v|, (2.D 

V V 

with a,j3 <G F2". Notice that in the self-dual basis these oper- 
ators factorize as 

Z« = of ® ■ ■ ■ <g> of" , Xp = a* 1 <g> • • • $ a b x " , (2.2) 

where a, = tr(a0,) and bj = tr(j30,-) are the corresponding ex- 
pansion coefficients for a and /3 in that basis. The single-qubit 
Pauli operators <7- and a T can be expressed in the standard ba- 
sis of the two-dimensional Hilbert space C 2 as 

ct z =|1}(1|-|0)(0|, cr r =|0)(l| + |l)(0|. (2.3) 

In addition, we have the commutation relation 

Z a X p = (-l)^X p Z a . (2.4) 

This is the discrete counterpart of the Heisenberg-Weyl group 
for continuous variables and the hallmark of noncommutativ- 
ity. Moreover, X a and Z a are related through the finite Fourier 
transform lfl9ll 

^=il(-l) tt(VV,) |v)(v'|, (2.5) 
V v ,v' 

so thatZa = &Za& S. 

We next recall fig] that the grid specifying the phase space 
for n qubits can be appropriately labeled by the discrete points 
(a, j3), which are precisely the indices of the operators Z a and 
Xp: a is the "horizontal" axis and j3 the "vertical" one. On 
this grid one can introduce a variety of geometrical structures 
with much the same properties as in the continuous case ll2"lll : 
the simplest are the straight lines passing through the origin 



(also called rays). These rays have a quite remarkable prop- 
erty: the monomials {ZaX^a} labeled by points of the same 
ray commute with each other, and thus have a common sys- 
tem of eigenvectors, which we shall label as |v,ju). Without 
going into details, they can be constructed as 

\v,p)=X v V ll \0)® — ®\0), (2.6) 

where is a finite rotation that changes the slope ji of the 
rays. Therefore, states with the same slope index fj, span an 
orthogonal basis. The explicit construction of can be found 
in Ref . |21J1 . In this way, both the state index v and the slope 
(i.e., basis) ji run over the 2" elements of F2L Of course, 
there is an extra basis of "infinite" slope (corresponding to the 
"vertical" axis) that cannot be obtained through a rotation: 

\v)=X v ^\0}®---®\0). (2.7) 

Let us enumerate for the time being these vectors by \v,k), 
where Ic runs the 2" values of in d2.61 l plus the extra basis in 
(I2.71 l. One can easily check that 

|<v,*|v' l *'>| 2 = 5,,,<w + ^(1 - So/) , (2.8) 

so they constitute a set of MUBs. In other words, the complete 
set of (2" + 1) mutually unbiased projectors 

^y,k = \V,k)(V,k\, (2.9) 

defines a complete scheme, in the sense that the measured 
probabilities 

p v>k = Tr(p^ Vjjt ) , £/> Vj * = 1 , (2-10) 

V 

determine completely the density matrix through 

2" + l 

P + l=II?v,i^. (2.H) 

k=l v 

Here, we have used Tr (with capital T) for the standard trace 
in Hilbert space. Note that the structure of this MUB set is 
preserved under any local unitary transformation, so any fac- 
torizable, complete basis can be chosen as a computational 
basis. 



III. OPTIMAL MEASUREMENT SCHEME UNDER 
PERMUTATIONAL SYMMETRY 

As heralded in the Introduction, we wish to design a mini- 
mal MUB tomographical scheme for systems that remain in- 
variant under all possible interchanges of its different parti- 
cles. This invariance could be simply stated as 

n pq pu pq = p, (3.i) 

where p 7^ q, p,q £ {1,...,«}. The elements Yl pq of the 
permutation group are known as swap operators, as they ex- 
change the states of the p-th and the g-th qubits; i. e., 

Tl pq \...,a p ,...,a q ,...) = \ ...,a q ,...,a p ,...) . (3.2) 
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It has recently been shown l22ll that any permutationally 
invariant n-qubit state defined via Eq. d3.lt can be written as 

n/2 

Ppi= PjPj R i- ( 3 - 3 > 

J Jmm 

The summation runs over different total spin numbers starting 
from y'nrin £ {0, 1 /2}, depending on whether n is even or odd. 
Pj is the density matrix of a j-spin state and pj are the associ- 
ated probabilities. In addition, Rj = 11 / dimJtfj and the factor 
dim Jdj comes from the degeneracy of the subspaces appear- 
ing in the decomposition of the total Hilbert space Jif = C 2 " 
in the form 

n/2 

J?= JtfjmXj. (3.4) 

J 7 min 

Here, Jfj is a spin Hilbert space of dimension dimJ^j =2j + 
1, while J£j are referred as multiplicative spaces that account 
for the different possibilities to obtain a spin j. One can show 
that JH 

^={n/2-j)-{n/2-j-l)- ^ 

This means that a permutationally invariant density operator 
only contains nontrivial parts of the spin Hilbert spaces, and 
there are no coherences between different spin states. In other 
words, any of these states can be parsed into a block-diagonal 
form that has been exploited in several contexts ll24ll . 

The crucial observation for what follows is that, at the level 
of field elements, the action of the permutation operator Yl pc/ 
on a state is tantamount of 

jU i->-jU+£tr(jU£), (3.6) 

where e = p + 9 q , with 9 p (0 q ) being the p-th (g-th) element 
of the self-dual basis. Since the field element addition is com- 
mutative, the operator is symmetric in p and q, as it should. In 
algebraic terms, we have then 

n M = £|jc+etr(ejf))(jc|. (3.7) 

K 

Using this field representation, we can check that the mutually 
unbiased projectors ( 12. 9\ transform as 

^[pqSP v,^n w = £P v+£tr(v£), P +£tr(v£) j (3-8) 

whence the density matrix in the tomographic representa- 
tion ( 12.1 \\ is transformed as 

Hpg (P + l)n w = ^/?V,/i^ v+£tr(V£).|i+£tr(V£) 

ii,v 

+ E/ 5 v^v+£tr(v£)- 0.9) 

V 

The last term is just the contribution from the basis (12.7) , 
which we split for notational simplicity. 



If we perform the change of variables 

v' = v + etr(ve), /t' = /t + £tr(;U£), (3.10) 
we can recast (13. 9\ in the form 

n p? (P + l)n M = Pv'+etr(v'e),/l'+etr(/l'e)^v',/l' 

+ E^v'+£tr(v'£)^v'- (3-11) 
v' 

Consequently, the invariance condition ( 13. U leads to the fol- 
lowing restrictions on the measured probabilities: 

Pv+etr(ve),n+etr(ne) = Pv,fi} Vfi. (3.12) 

Obviously, the probabilities p v )1 should be also invariant un- 
der all consecutive index permutations. 

IV. PHYSICAL DISCUSSION 

The above basic expression can be given a transparent phys- 
ical meaning. Indeed, let us expand v and fi in the self-dual 

basis 

n n 

M = £m i -0 i , V = (4.1) 

1=0 i=0 

with nij , rti £ Z2 and analogous expansions for the transformed 
indices in (13.10b . For a given £ = 9 p + q , one can check that 
m'j = mu except for m' p = m q and m' q = m p , and similarly for 
n'j. That is, a change of the index v of the states in a basis 
simply results in a reshuffling of its states. Therefore, such 
transformations do not give any new tomographic projectors 
for a permutationally invariant state. 

The transformation of fi implies that measurements by 
MUBs corresponding to fi indices with the same number of 
non-zero components in the self-dual basis [the length of the 
word I jit I corresponding to the binary string (ao,a\ r ..,a„)] 
give the same information. In short, the projectors labeled by, 
e.g., jU = (1, 1,0, . . . ,0) and \l' = (0, . . . ,0, 1, 1) are equivalent. 
The computational basis, associated with fX=0, automatically 
satisfies ( 13.121 i for all v. Similarly, the X basis also satis- 
fies (13.12b because it has no fi dependence ([n w , J^] = 0). 
Therefore, these two bases remain invariant under any qubit 
permutation. This allows us to count the total number of mea- 
surements needed for a complete reconstruction of the density 
matrix, which is just n + 2. This result could be expected, 
for the Hilbert space dimension of the permutation invariant 
system is n+ 1. 

Since the permutation group acts simultaneously on both 
indices fi and v, there are different orbits of equivalent prob- 
abilities that are defined not only by the length |ju| but also 
by the mutual symmetry properties of the indices representing 
the number of the basis and the element in each basis. In par- 
ticular, for jx ^ each orbit representative is labeled by three 
lengths m=|/i|,/ = |v|,5 = |^i + v|, i.e. p v/1 =p(m,l,s). For 
the computational and the Fourier bases the orbits are charac- 
terized only by |v|; for instance, p v .o = po(l)- Accordingly, 
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in each basis not all the probabilities should be measured, but 
only those that belong to different orbits, which leads to a re- 
duction of the experimental errors. Since for given m and /, s 
runs from \m — 1\ to min(m + l,2n — m — l,n) in steps of two, 
the number of orbits turns out to be 1 +n(n 2 +6n + 17)/6. 
Bearing in mind the normalization condition (12.10b . we find 
that there are n(n 2 + 6n + 17) / 6 independent probabilities 
p(m,l,s), which completely define the density matrices ap- 
pearing in the decomposition (13.3b . Projectors corresponding 
to the same probabilities are given by the condition ( 13 . 1 2b . 
The final reconstruction takes the form 

n n min(m+/,n) 

p+i = ££ E p( ffl .'. f )EE^(/, S )^H 

m=0/=l s =\m-l M V 

+ 1 Mi) L ^vco,o + 1 p(i) L + Lp°^ ^ 

1=1 v i=o v fiv 

(4.2) 



where the sum on jx and v run over all the field elements such 
that \fi(m,s)\ =m,l= |v(/,j)|, s — |ju + v|, and po^ = 1 - 
EvPv,ii, and p = 1 -£ v p v . 



For instance, for two qubits, the field F 2 2 has the primitive 
element defined by the irreducible polynomial 9 2 + + 1 =0. 
Therefore, Q x = 6 and 2 = 2 , so that 3 = dd 2 = 0(1 + 
0) = 0\ + 02- In this case, only measurements in the bases 
with pi = 0i (or ji — 02) and fi = 0\ + 02 (apart from measure- 
ments in the computational and X bases) are required. The 9 
independent measured probabilities [pe,,o and pg, + e- Q, from 
theZbasis, pe u ei, and pe 1+ e,,ei from basis 1, /?e,,e 1+ e 2 

and pe l +e 2 .e [ +e 2 from basis 3, and pq { and pe i +e 1 from the X 
basis] are representatives of the equivalent probabilities orbits. 
This selection gives an explicit reconstruction form that reads 
as 



p + 1 = p 8li o (^Vo + ^*e 2 ,o) +P0i +e 2 fi3 B e l +e 2 ,o 

+ Pe 1 ,e 1 {^e u Bi + ^e 2 ,e 2 ) + pe 2 .e l (^Ve, + ^ , e ll 2 ) +/>e 1 +e 2 ,e I (^ei+e^e, +^81+82,%) 
+ Pe l ,e l +e 2 (^ 2 e 1 ,e 1 +e 2 + ^fl^fli+fli) + Pe l +e 2 ,e l +e 2 ^'e l +e 2 ,e l +e 2 
+ Pe l (^e, + ^e 2 ) +Pe 1 +e 2 ^ > e 1 +e 2 

+ Po,o^ a o,o + /'o,e l ^ l o.ei + /?o,e 2 ^o,e 2 +Po,8i +02^0,0! +e 2 +P0&0, (4.3) 



r 



where po./j = 1 — Lv Pv,n an d thus can be derived from the 
nine independent, measured probabilities. Similarly, pQ = 1 — 
2 P0i -pe 1 +e 2 - 

For this problem the computational basis is 



{X v |0)} = {|v)> 



(4.4) 

The three remaining bases (apart from a normalization factor) 
are 




{Xv|6i)} 



{X v |0 2 >} 



{X v |0 1 + 2 )} = 




while the one corresponding to (12.7b turns out to be 
1 

|v)} = 




(4.6) 



Of these five MUBs, only, e.g., {|v)}, {X v |0i)}, {*v|0i + 
62)}, and the {|v}} are needed to tomographic ally recon- 
struct a permutationally invariant two-qubit state. If we per- 
mute the second and the third qubit (and the state to be to- 
mographed would not change due to such permutation) it is 
readily seen that the permuted basis {X v |02)} becomes the 
non-permuted basis {X v 1 9\ ) } (but with the middle two vectors 
interchanged). Hence, the two bases extract identical informa- 
tion from the state, and hence one of them can be disregarded. 
They both have one nonzero component in the self dual basis 



TABLE I. Allowed values of m, 1 and s for the 24 independent orbits 
in the three-qubit case. The tilde indicates that the corresponding 
probabilities are measured in the X basis. The last row (denoted #) 
gives the number of (equivalent) probabilities in each orbit. 

m OOOOOO661 1 1 1 1 12222223333 

/ 012301230112230112230123 

s 012301231021322130313210 

# 133113313366333633631331 
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and are therefore directly related by a permutation as shown 
by (EE). 

Before we conclude, let us briefly address the case of three 
qubits. In Table Q] we give the values of l,m,s for the 24 inde- 
pendent orbits (all in all, we get 72 probabilities). Taking into 
account that 5 probabilities (each one defining an orbit) can be 
determined from the normalization condition (12.10) [for ex- 
ample, we can fix p(0,0,0), p(0,0,0,), p(l,0,l), p(2,0,2), 
and p(3,0,3)], we arrive at 19 orbits that determine any sym- 
metric density matrix. 

V. CONCLUSIONS 

We have developed a method to generate a minimal set of 
MUBs needed to tomographically reconstruct a state consist- 
ing of n qubits, when the state is invariant under the permu- 
tation of the qubits. Such a state spans an n + 1 dimensional 
Hilbert space. Consequently the smallest set of bases one can 
hope to use is n + 2, and indeed our method provides a mini- 
mal set. 



MUBs are not strictly necessary to reconstruct such a state, 
but they have the advantage of capturing maximally different 
aspects of the state. Moreover, as the bases constitute com- 
plete sets of states in the 2" -dimensional space of n qubits, 
they can in principle be implemented as von Neumann mea- 
surements and not as individual projectors or positive operator 
valued measures. The price is that the MUB projectors are for 
the most part highly entangled, so their experimental imple- 
mentation can be difficult. 
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